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We report valence and conduction band alignments and offsets for heterojunctions between 
CdCr2Se4, an n-type ferromagnetic semiconductor, and the non-magnetic materials Si and GaAs, 
evaluated using density functional theory. We explore numerically the impact of different interface 
features on the type of band alignment and the magnitude of the offsets. For example, we find it is 
energetically favorable to deplete Cr atoms from the layers at the interface; this also leads to band 
alignments smaller in magnitude compared to those obtained for Cr-rich interfaces and ideal for 
electrical spin-injection into either Si or Ga-terminated GaAs substrates. 



I. INTRODUCTION 

Spin "injection" refers to the induced transport of 
spin-polarized electrons or holes from one material into 
another. To date, two basic approaches to injection 
have been explored: photoinduced generation of spin- 
polarized carriers in one material, which thpp diffuse 
across an interface into a second material ;ETbil and di- 
rect electrical injection.. via application of a bias volt- 
age across an interfaceJjLl In both approaches, the het- 
erojunction band offset represents an important quan- 
tity which may place fundamental constraints on the ex- 
pected efficiency of the injection process. In this work 
we report a first-principles numerical study of band off- 
sets for realistic heterojunctions between magnetic and 
nonmagnetic semiconductors. 

We will refer to the magnetic and nonmagnetic sides 
of the heterojunction as the "source" and "sink," respec- 
tively. Our choice of the source and sink materials are 
guided by a few simple principles. First, since spin life- 
times are generally longer for electronsa and since elec- 
trons generally have higher jpobilitics and smaller spin- 
orbit interaction than holes ,□ we focus on source materi- 
als which can supply spin-polarized electrons rather than 
holes. Second, we restrict our attention to semiconduct- 
ing source materials so as to avoid the well-known "con- 
ductivity mismatch" problem that plagues direct injec- 
tion from a metal into a semiconductor £3 Third, because 
lattice mismatch leads to dislocations (which act as spin 
scatterers) , we consider sink materials with a reasonably 
close lattice match to the source material. A particularly 
interesting candidate source material is the chromium 
chalcogenide spinel CdCr2Se4, an n-type ferromagnetic 
semiconcLuctor with a relatively high Curie temperature 
of 130 K.EJ For sink materials, we consider Si and GaAs, 
both of which can be made n-type with appropriate dop- 
ing and are well lattice matched to CdCr 2 Se4. 

In this work we use density functional theory to cal- 
culate band alignments and offsets between CdCr 2 Se4 
and Si or GaAs. We take the view that although the 
band alignment may ultimately place constraints on the 
injection efficiency, one may be able to "tailor" the inter- 
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FIG. 1. Various types of band alignments in a semicon- 
ducting heterojunction. Valence and conduction band edges 
are denoted vb and cb, respectively. 

face during the growth process, thereby improving align- 
ments that are otherwise unfavorable for injection. It 
is this variability in the alignments and offsets, with re- 
spect to various features of the interface, that we explore 
numerically in this work. Of the three possible types of 
alignment (see Fig. [I]), the two most favorable for spin- 
injection of electrons are Type I or IIA. Less favorableis 
Type IIB alignment, which leads to Zener breakdownlla 
(valence electrons of the source material tunnel into the 
conduction bands of the sink). Type IIB alignment also 
gives rise, through band bending, to confinement of holes 
in the source material in the layers adjacent to the inter- 
face, which would then act as strong spin-scattering cen- 
ters and presumably degrade spin transport across the 
interface. 



II. THEORY OF BAND OFFSETS 

In this work we evaluate the band offsets usras-jneth- 
ods pioneered by Van de Walle and Martin,t2rtZl here 
extended to include spin splitting. For completeness we 
review the essentials of that approach. We take the sink 
material to be on the left (L) side of the heterostruc- 
ture, and the source material to be on the right (R). The 
valence band offset (VBO) is the relative alignment of 
the L and R valence band maxima (VBM) far from the 
interface, and is conveniently written as 
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The various contributions to the valence band 
offset, -Evbo- Vavg is the average value of the effective poten- 
tial in L and R far from the interface, to which the valence 
band maxima are referenced. Also shown are the level shifts 
arising from strain in the overlayer (R), as well as spin-orbit 
(SO) splitting in L and R. Only one spin channel is depicted. 
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Here Ai5band and AV avg are generally referred to as the 
band lineup and potential lineup, respectively. Ai?band 
is determined entirely by the electronic structure of the 
constituent bulk materials, while AV^ V g accounts for the 
effects of the interface between those materials. For sim- 
plicity we do not show spin indices in Eq. |l|; each spin 
channel is treated independently. 

Although the band lineup, Ai^band, is a purely bulk 
quantity, its evaluation in general depends on the na- 
ture of the interface. For example, lattice mismatch be- 
tween the substrate and overlayer will introduce tetrag- 
onal strain into an overlayer that has a nominally cubic 
lattice. This strain will in general lead to level shifts 
and band splittings that must be correctly included when 
computing the bulk band lineup. For the purposes of this 
work we consider the nonmagnetic sink materials to be 
unstrained substrates on which overlayers of CdCr2Se4 
are grown epitaxially. Thus the overlayers experience an 
in-plane tensile strain (1.3% and 3.3% for Si and GaAs, 
respectively) that results in a contraction of the spacing 
between the overlayers. To fully account for these strain 
effects we consider bulk unit cells for CdCr2Se4 with a 
tetragonal distortion a\ = a-i = 2a su bstratc ^ ^3- 

We also include fully the effects of the spin-orbit (SO) 
interaction on the band structure, by performing fully rel- 
ativistic calculations for the bulk phases. With strain and 
spin-orbit effects thus included, the band lineup takes the 
form 
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Here, -Evbm is the position of the valence-band maxi- 
mum, and V avg is the cell-averaged Kohn-Sham potential 
(which provides a convenient reference energy). AE^o 



is the shift in the VBM associated with spin-orbit in- 
teraction. Figure ^ shows an an energy level-diagram 
including all the various terms defined above. 

The second term in Eq. [j], the potential lineup AVa Vg , 
is determined by the electronic structure of a specific in- 
terface. We use periodic supercells to represent the het- 
erojunction, with supercell periods sufficiently large to 
ensure that the interfaces do not interact. The poten- 
tial lineup is the difference in the locally averaged po- 
tential on the two sides of the interface, far from the 
interface itself. To compute these local averages we use 
the slab-averaging technique discussed in Rcf. [r| One 
first computes the planar average, V(z), of the total po- 
tential, V(r), by averaging across planes parallel to the 
interface. V(z) contains atomic-scale variations whose 
local periodicity on the two sides of the interface may in 
general may be different — even far from the interface — 
because of lattice mismatch. The correct local averages 
can be obtained by performing a "slab average" over a 
sliding window. This slab averaging is applied-twice, in 
order to allow for different L and R windows £j£j 
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Here Xl and A_r denote the periodic repeat distances in 
the direction perpendicular to the interface for the L and 
R materials . Far from the interface, the slab-averaged 

potential, V(z), is constant. The potential lineup is then 
given by 



AV avg = V(z L ) - V(z R ), 



(4) 



where z^ and zr denote the midpoints of the L and R 
slabs. Finally, having computed the VBO directly from 
Eq. |l|, the conduction band offset (CBO) is simply 
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where -E g a P is the band gap for L (R). For reasons dis- 
cussed below (see Sec. ^), we evaluate Eq. || using exper- 
imental band gaps rather than our theoretical values. 



III. INTERFACE MODELS 

The crystal structure of bulk CdCr 2 Se4 is shown in 
Fig. 0. In this work we consider only (001) interfaces, 
for which tetrahedral bonding across the interface can be 
maintained in a particularly simple manner. In Fig. |^, 
different (001) basal planes can formed by truncating the 
spinel structure at different atomic planes. Only two 
chemically inequivalcnt planes can be so formed: one 
containing only Cd atoms, and one containing Se and 
Cr atoms in a 2:1 ratio. Such a representative plane con- 
taining Se and Cr atoms is indicated in Fig. ||. Within 
these planes, the atomic density in the Cd plane is only 
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FIG. 3. Crystal structure of CdCr2Se4 in the bulk spinel 
phase. Representative Cd, Cr, and Se atoms are labeled. The 
Se atoms form the vertices of octahedra (some of which are 
shown) centered at each Cr site. A (001) plane containing Cr 
and Se atoms is also indicated. 

one fourth that of Si or GaAs (ignoring the lattice mis- 
match), which suggests that a 2x2 reconstruction of the 
terminating Si layer would result for interfaces with the 
Cd plane. On the other hand, for interfaces with the 
mixed (Se,Cr) plane a simple unreconstructed interface 
is quite plausible, and so we limit ourselves to variants of 
this basic interface geometry. 

For interfaces with Si, the atomic alignment at the in- 
terface is shown schematically in Fig. |^. For the case of 
Si we consider only atomically abrupt interfaces. At the 
interface boundary, every Si atom is tetrahedrally coor- 
dinated (to two Si and two Se atoms). The role of Cr 
atoms at the interface is not immediately obvious from 
these preliminary considerations, but will be analyzed in 
detail in Sect. |v| 

For interfaces with GaAs we consider both abrupt and 
partially intermixed interfaces. In both cases we com- 
pute offsets for both Ga-terminated and As-terminated 
substrates. Fig. || shows an abrupt As-terminated sub- 
strate; Ga-terminated substrates are analogous. 

Although abrupt interfaces are conceptually attractive, 
for a polar material like GaAs they may not be physi- 
cally plausible. In particular, for interfaces between po- 
lar and non-polar materials, abrupt termination leads ip* 
evitably to the formation of long-range dipole fields.tS 
Such a dipole field is energetically unfavorable over large 
distances, and will typically lead to an atomic rearrange- 
ment near the interface in order to reduce such effects. 
For polar/nonpolar interfaces, a particularly simple way 
of eliminating the dipole field is to intermix the polar and 
nonpolar constituents at the interfacial layer in order to 




FIG. 4. Schematic of an abrupt [001] interface between 
CdCr2Se4 and Si. Representative Cd, Cr, Se and Si atoms 
are denoted, as are the outlines of the octahedron defined 
by the Se atoms. Some aspects of this figure (bond lengths, 
angles, etc.) are not drawn to scale. Note that not all atoms 
in this figure lie in the same plane. 




FIG. 5. Schematic of an abrupt [001] interface between 
CdCr2Se4 and GaAs. 

maintain local charge neutralityQ We have followed this 
approach for interfaces with GaAs, by appropriately sub- 
stituting Ga or As atoms for Se atoms at the interface in 
such a way as to maintain local neutrality. Fig. |^ shows 
an example of this substitution for an As-terminated in- 
terface. 



IV. COMPUTATIONAL DETAILS 

For both the interface and bulk calculations we use 
the plane-wave pseudopoteniial method with norm- 
conserving pgeudopotentials,E3 as implemented in the 
ABINIT code.EJ In both methods we use the local-spin- 
density approximation (LSDA) in the form of Perdew 
and WangnJ We have carefully checked the accuracy of 
the pseudopotentials by comparing to results for bulk 
CdCi'2Se4 obtained with the linearized augmented planc- 
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FIG. 6. Schematic of partially intermixed [001] inter- 
face between CdCr2Se4 and GaAs. From simple counting 
arguments As-Se bonds contribute an extra | electron while 
Ga-Cd bonds contribute an excess § hole. Thus substitution 
of Ga for Se in the interface layer of plane (a) [cf. Fig. ^|(a)] 
maintains local charge neutrality over the first few interface 
layers. 

TABLE I. Structural and magnetic properties of bulk 
CdCr2Se4, as predicted by the pseudopotential plane- wave 
and LAPW methods and as measured experimentally. The 
four columns are lattice constant, bulk modulus, magnetic 
moment, and normal- vs. inverse-spinel total energy. 



Method 


a (A) 


B (GPa) 


M (/x B ) 


AE S (eV) 


PSPW 


10.68 


96.2 


6.0 


-2.72 


LAPW 


10.59 


93.9 


6.0 


-1.97 


Exp. 


10.72 a 


76.9 b , 43.4 C 


5.4 b , 5.6 d 





a Ref. 
b Ref. 
c Ref. 
d Ref. 



wavepXLAPW) method as implemented in the wien97 
codeB Table m summarizes the predictions of both meth- 
ods for equilibrium lattice constant, bulk modulus, and 
magnetic moment, along with the available experimental 
values. The pseudopotential and LAPW lattice constants 
agree to better than 1%, the bulk moduli to better than 
3%, and the magnetic moment essentially exactly. The 
agreement with experimental values is equally good with 
the exception of the bulk modulus (for which the spread 
in experimental values is very large) . We have also com- 
puted the difference in total energy between the normal- 
and inverse-spinel forms of CdCr2Se4, 

A_Z? S ^normal -^inverse j (6) 

and find that the spinel structure is energetically pre- 
ferred (see Table |) within both PSPW and LAPW meth- 
ods. Thus for the rest of this work we consider only the 
normal spinel phase. 

The interface results we report here were obtained with 
supercells containing five atomic layers of Si or GaAs 



and five layers of CdCr2Se4. These thicknesses are suf- 
ficient to ensure that the interfaces do not interact, as 
confirmed by the convergence (discussed below) of the 
planar-averaged potentials to their bulk values within the 
first few layers of the interface. The coordinates of all the 
atoms were completely relaxed until all force components 
were smaller than 10~ 3 Ha/bohr. By simultaneously op- 
timizing the supercell dimension normal to the interface, 
we impose no artificial constraints on the interlayer spac- 
ing throughout the supercell. As described above, how- 
ever, we do constrain the in-plane lattice constant to that 
of Si or GaAs. The resulting tensile strain experienced 
by the CdCr 2 Se4 overlayers gives rise to a decrease in 
the spacing between bulk-like CdCr2Se4 layers, consis- 
tent with total energy calculations for the bulk material 
with a tetragonal distortion. The cumulative effect of all 
these physical relaxations (relative to the ideal geometry) 
on the absolute band offsets is typically of the order of 
0.1-0.3 eV, and so must be fully included for an accurate 
description of the alignments. 

For all the pseudopotential calculations, a plane-wave 
cutoff of 70 Ry was used, which is sufficient to converge 
all band states relevant to the offset-calculations to within 
100 meV. A 4x4 Monkhorst-PackcH sampling of the sur- 
face Brilloun zone was used for the supercell calculations, 
equivalent to 8x8 sampling of the primitive Si or GaAs 
surface zone. 



V. BULK CdCr 2 Se 4 BAND STRUCTURE 

The electronic structure of bulk CdCr2Se4 has been 
the subject of a number of experimental studies. Pho- 
toluminescence experiments give the fundamental gap — 
the energy difference between the highest p-like valence 
state and lowest s-like conduction state — to be 1.8 eV 
in the ferromagnetic ph.ase.IEl In contrast, optical studies 
reveal an absorption edge only 1.3 eV above the valence- 
band edge, suggesting the existence of d-like conduction 
states in the fundamental gap.ES This view is supported 
by transport measurements showing relatively low mobil- 
ity, consistent with conduction taking place within nar- 
row d-like bands .o 

Our results for the band structure of bulk CdCr2Se4 
are show in Fig. [?]. We find that the ground state 
is indeed a ferromagnetic semiconductor, with a mag- 
netic moment of 6.0 hb per formula unit. Our results 
are in very rgpod agreement with other recent LAPW 
calculations," but are slightly different from results 
obtained with the linearized-muffin-tin-orbital (LMTO) 
method, which gives semimctallic band structure in both 
spin channels-and a slightly higher moment of 6.2 /ie per 
formula unit 13 

Figure fj] shows the density of states for bulk CdCr2Se4 
decomposed into angular character and projected onto 
atomic sites. Each Cr site contributes a moment of 
3/is to the total moment, corresponding to four spin-up 
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FIG. 7. Electronic structure of bulk CdCr2Se4 in the fer- 
romagnetic phase for the majority (upper panels) and mi- 
nority (lower panels) spin channels. The arrows labeled "fg" 
marks the s-p fundamental gap. The valence band maximum 
is taken to be the energy zero. 



and one spin-down electron. Our predicted fundamental 
gaps, indicated in Fig. [?| by arrows, are found to be 2.42 
and 3.05 eV for the majority and minority spin channels 
respectively; these results are in very good agreement 



with those reported in Ref. 34 



For spin injection, the relevant energy gap is the band 
gap between valence and conductions states. Our re- 
sults predict the band gap to be very small, of order 
0.3 eV in both spin channels; this value is much smaller 
than the experimental optical absorption edge at 1.3 eV. 
Since there is experimental evidence suggesting that con- 
duction takes place_in_the unoccupied d levels probed by 
optical absorption,E§E3 we simply adopt the experimental 
value of 1.3 eV for the band gap of CdCr2Se4 appearing 
in Eq. For similar reasons, we use experimental values 
of the band gaps of Si and GaAs.LJ 



VI. CALCULATED BAND OFFSETS 

We begin this section by illustrating the slab-averaging 
procedure used to determine the potential lineup across 
an interface. Figure || shows the planar-averaged and 
slab-averaged potentials for the case of CdCr 2 Se4/Si. 
The five layers of Si are in the middle portion of the 
figure, the atomic planes corresponding to minima in 
the planar-averaged potential; the three maxima in the 
CdCi'2Se4 portion of the cell correspond to Cd planes. 
The slab-averaged potential is quite flat at the mid- 
points of the Si and CdCr 2 Se4 slabs, demonstrating that 
these slabs are sufficiently thick to converge the potential- 
lineup contribution to the band offset. Figure || also 
shows the planar-averaged potential for bulk Si. The 
rapid convergence of the planar-averaged potential in the 
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Cd Si Si Si Si Si Cd Cd 

Planar (dashed curve) and slab-averaged (solid 
curve) potentials in a CdCr2Se4/Si heterojunction for the ma- 
jority spin channel (minority spin channel is similar). The po- 
tential lineup is also indicated. The planar-averaged potential 
for bulk Si is shown for reference (dotted curve). 



supercell to its bulk value is readily apparent, and is sim- 
ilar to that found in previous studiesJ13 

Our calculated valence- and conduction-band offsets 
for the five interfaces described in Sec. [II are summa- 



rized graphically in Fig. ^ and tabulated numerically in 
Table ^ In Table Q we show only offsets for the valence 
bands since only these were evaluated completely from 
first principles. For all five interfaces studied, we find 
that the CdCr2Se4 conduction-band edges for both spin 
channels sit well above the conduction-band edge of the 
substrate — independent of the substrate material (Si or 
GaAs), termination layer (Ga or As), and the atomic in- 
termixing described in Sec. III. The band alignments are 
all Type IIB for the majority spin channel, and either IIA 
or weakly IIB for the minority channel. Thus, according 
to the criteria outlined in Sec. |, these band alignments 
are not favorable for efficient spin injection. 

A major goal of this work is to identify specific in- 
terface features that strongly affect band offsets, and to 
evaluate whether such features can be externally con- 
trolled, for example, during the growth process. For sev- 
eral of the interfaces we have studied, we have identified 
one such feature: the concentration of Cr atoms at the 
interface. Our motivation for studying the effects of vary- 
ing the interfacial Cr content comes from analyzing the 
environmental dependence of Cr removal energies. For 
example, the calculated formation energy for a Cr va- 
cancy in bulk CdCr 2 Se4 is 3.27 eV, but for Cr atoms in 
the interface layer of a CdCr2Se4/Si interface the forma- 
tion of a single vacancy is actually exothermic, with an 
energy gain of 1.60 eV/vacancy£2 

To explore the impact of interfacial Cr content on band 
offsets, we consider the limiting case of total depletion of 
Cr from the mixed (Se,Cr) layer at the interface. Since 
the Cr atoms do not participate directly in the tetrahc- 
dral bonding across the interface boundary, their removal 
leads to relatively minor rearrangements of the interface 
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TABLE II. Valence band offsets for hetero junctions of Si 
and GaAs with CdCr2Se4. Offsets are in eV. The notation 
A|(B,C) denotes a substrate terminated by atoms of chemical 
type A adjacent to interface layers containing atoms of type 
B and C. For example, Fig. ^ is denoted Si|(Se,Cr), whereas 
Fig. | is denoted As|(Se,Cr,Ga). 



Heterojunction 


VBO T 


VBOi 


Si|(Se,Cr) 


1.38 


0.79 


Si (Se) 


0.37 


-0.14 


As|(Se,Ga,Cr) 


1.87 


1.29 


As (Se,Ga) 


1.78 


1.21 


As (Se,Cr) 


2.09 


1.55 


As (Se) 


1.99 


1.47 


Ga|(Se,As,Cr) 


1.87 


1.32 


Ga (Se,As) 


1.30 


0.78 


Ga (Se,Cr) 


1.95 


1.41 


Ga (Se) 


1.46 


0.94 



structure — which we nevertheless fully account for by re- 
optimizing the geometry. The energetics of removing this 
Cr layer vary somewhat from interface to interface, but 
even in the limiting case of a Cr-rich growth environment, 
the total depletion of Cr from the final interface layer 
is thermodynamically favorable for all five interfaces, as 
shown in Table 



III 



The band alignments and offsets for Cr-depleted inter- 
faces are summarized in Fig. 
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Compared to the non- 
depleted interfaces, all of the offsets are reduced. The 
reduction is remarkably large for the interface with Si 
(about 1 eV) , smaller for both abrupt and intermixed Ga- 
terminated GaAs interfaces (about 0.5 eV), and smallest 
for As-terminated GaAs interfaces (about 0.1 eV). In two 
cases (Si, and Ga-terminated GaAs) the reduction is ac- 
tually sufficient to change the alignment type from IIB 
to IIA. 

Finally, we emphasize that the question of whether the 
concentration of interfacial Cr can be externally con- 
trolled remains open. Although the energetics clearly 
favors Cr-depleted interfaces, kinetic barriers may be 
more important — and perhaps decisive — at the relatively 
low growth temperatures typical of molecular-beam epi- 
taxy. Theoretical studies of complex interfaces are no- 
toriously difficult, not the least because of uncertain- 
ties in the interface morphology. The present work 
clearly demonstrates the possible dangers of assuming 
an "idealized" interface structure, and should make clear 
that an understanding of the systematics of band offsets 
requires — at the very least — detailed microscopic calcu- 
lations. TABLE 



TABLE III. Energy gained-by depleting Cr from the inter- 
face layer, in eV per vacancy." The labels refer to the parent 
heterojunction before depletion of Cr. 



Heterojunction 



AE 



Si|(Se,Cr) 
As|(Se,Ga,Cr) 
As|(Se,Cr) 
Ga|(Se,As,Cr) 
Ga|(Se,Cr) 



0.69 
0.81 
0.38 
0.12 
0.91 



Si | (Se,Cr) As | (Se,Ga,Cr) Ga | (Se,As,Cr) As | (Se,Cr) Ga | (Se,Cr) 




FIG. 9. Valence and conduction band offsets for hetero- 
junctions with Si and GaAs for majority and minority spin 
channels. Spin-orbit splitting of the valence bands is denoted. 
The notation labeling the interfaces is described in the cap- 
tion for Table II. 



VII. SUMMARY 

In summary we have applied first principles methods 
to determine the absolute valence band offsets for het- 
erojunctions between CdCr2Se4 an n-type ferromagnetic 
semiconductor, and Si or GaAs. By using experimen- 
tally determined band gaps, we have also determined the 
type of band alignments in these heterojunctions, and 
find favorable alignment for electrical injection for Cr- 
depleted heterojunctions with Si and with Ga-terminated 
GaAs substrates. These findings suggest that CdCr 2 Se 4 
may be a promising candidate for the injection of spin- 
polarized electrons. Finally, we point out that although 
the type Type IIB alignments we have found are not 
favorable for spin injection, such alignments are ideal 
for spin-polarized resonant interband tunneling diodes in 
which CdCr2Se4 would function as a spin-selective quan- 
tum well. 



VIII. ACKNOWLEDGEMENTS 



One of the authors (J. M.S.) acknowledges the Na- 
tional Research Council for support during this work in 
the form of a postdoctoral fellowship. We thank An- 
dre Petukhov, Steve Hellberg, Berry Jonker and Aubrey 



G 



Si | (Se) As | (Se,Ga) Ga | (Se,As) As | (Se) Ga | (Se) 




(1991). 

21 The ABINIT code is a common project of the Universite 
Catholique de Louvai n, Corning Incorporated, and othe r 



vb )IH vb 

FIG. 10. Valence and conduction band offsets for hetero- 
junctions with Si and GaAs, in which the interface layers has 
been depleted of Cr. 



Hanbicki for many useful conversations during the course 
of this work. This work was funded in part by DARPA 
and ONR. Computational work was supported in part 
by a grant of HPC time from the DoD Major Shared 
Resource Center ASCWP. 



1 S. Koshihara et al, Phys. Rev. Lett. 78, 4617 (1997). 

2 I. Malajovich et al, Phys. Rev. Lett. 84, 1015 (2000). 

3 I. Malajovich et al, Nature 411, 770 (2001). 

4 R. Fiederling et al, Nature 402, 787 (1999). 

5 Y. Ohno et al, Nature 402, 790 (1999). 

6 B. T. Jonker, et al, Phys. Rev. 62, 8180 (2000). 

7 A. T. Filip et al, Phys. Rev. B 62, 9996 (2000). 

8 J. M. Kikkawa and D. D. Awschalom, Phys. Today 52, 33 
(1999). 

9 W. A. Harrison, Electronic Structure and the Properties 
of Solids (W. H. Freeman and Company, San Francisco, 
1980), chap. 6. 

10 G. Schmidt, D. Ferrand, L. W. Molenkamp, A.T. Filip, and 
B. J. van Wees, Phys. Rev. B 62, R4790 (2000). 

11 H. W. Lehmann, Phys. Rev. 163, 488 (1967). 

12 See, for example, Principles of Electrical Engineering Ma- 
terials and Devices (McGraw-Hill, Boston, 2000), chap. 6. 

13 C. G. Van de Walle and R. M. Martin, Phys. Rev. B 34, 
5621 (1986). 

14 C. G. Van de Walle and R. M. Martin, Phys. Rev. B 35, 
8154 (1987). 

15 C. G. Van de Walle, Phys. Rev. B 39, 1871 (1989). 

16 A. Franciosi and C. G. Van de Walle, Surf. Sci. Rep. 25, 1 
(1996). 

17 A. Qteish and R. J. Needs, Phys. Rev. B 45, 1317 (1992). 

18 H. M. Al-Allak and S. J. Clark, Phys. Rev. B 63, 033311 
(2001). 

19 W. A. Harrison et al, Phys. Rev. B 18, 4402 (1978). 

20 N. Troullier and J. L. Martins, Phys. Rev. B 43, 1993 



contributors (URL http://www.pcpm.ucl.ac.be/abinit). 
This codep^elies on an efficient Fast Fourier Transform 
algorithm,E3 an ao\ajptation of the band-by-band conjugate- 
gradient method,H-i and a potential-based conjugate- 
gradient algorithm .cfl 

22 S. Goedecker, SIAM Journal on Scientific Computing 18, 
1605 (1997). 

23 M. C. Payne et al, Rev. Mod. Phys. 64, 1045 (1992). 

24 X. Gonze, Phys. Rev. B 54, 4383 (1996). 

25 J. P. Perdew and Y. Wang Phys. Rev. B 45, 13244-13249 
(1992). 

26 P. Blaha, K. Schwarz, P. Dufek, and R. Augustyn, 
WIEN95, Technical University of Vienna, 1995 (improved 
and updated UNIX version of the original copyrighted 
WIEN code, P. Blaha, K. Schwarz, P. Soratin, and S. B. 
Trickey, Comp. Phys. Commun. 59, 399 (1990)). 

27 P. K. Baltzer, H. W. Lehmann, and M. Robbins, Phys. 
Rev. Lett. 13, 493 (1965). 

28 B. Batlogg, M. Zvara and P. Wachter, Solid State Commun. 
28, 567 (1978). 

29 A. P. Galdikas et al, Sov. Phys. Solid State 31, 1228 
(1989). 

30 P. K. Baltzer et al, Phys. Rev. 151, 367 (1966). 

31 S. S. Yao et al, Phys. Rev. Lett. 46, 558 (1981). 

32 M. D. Coutinho-Filo and I. Balberg, J. Appl. Phys. 50, 
1920 (1979). 

33 N. Shanthi, P. Mahadevan and D. D. Sarma, Journ. Sol. 
State Chem. 155, 198 (2000). We note that those authors 
correctly point out a clear signature of a gap in their pre- 
dicted DOS (see their Fig. 8); however, the corresponding 
band structure shows this gap to be extremely small or 
zero. 

34 A. Continenza et al, Phys. Rev. B 49, 2503 (1994). 

35 P. Lederer, M. Heritier, and F. Brouers, Solid State Com- 
mun. 11, 1759 (1972). 

36 Numerical Data and Functional Relationships in Sci- 
ence and Technology, edited by K.-H. Hellwege, Landolt- 
Brnstein, New Series, Group III, Vol 17, edited by O. 
Madelung, M Schulz, and H. Weiss (Springer, Berlin, 1982). 

37 H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 
(1976). 

38 To evaluate these formation energies we have used the en- 
ergy per atom of bulk Cr as the chemical potential of Cr. 
Although the conditions of the growth process may lead to 
a different numerical value the bulk energy per atom is the 
thermodynamic upper bound, and hence any lower value 
leads to a larger energy gain for Cr depletion. 



7 



